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■ To clarify the origin of anomalous behaviors in Pb/S^RuCU/Pb junctions in terms of the pairing 

symmetry, a theory of the d.c. Josephson current in s-wave superconductor / p-wave superconductor 
/ s-wave superconductor junctions is developed. Calculated results on the temperature dependence 
of the critical Josephson current exhibit non-monotonous behaviors when the thickness of the p-wave 
superconductor is comparable to the coherence length. The consistency between present results 
with recent experimental measurement supports the possibility of a unitary p-wave pairing state in 
Sr 2 Ru0 4 . 
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Recently, physical properties on Sr2Ru04 have been studied vigorously both in the superconducting states and 
^ hi the normal state Jl]]. There have been presented several experimental results as well as theories that support an 
unconventional superconducting states i.e. triplet p-wave pair potential, in superconducting Sr2Ru04 |^-^|. In order 
to identify definitively thepairing symmetry, phase sensitive measurements utilizing the Josephson junctions give the 
most useful information JTyJll[] . Standing on this view point, it is an interesting problem to develop a theory of the 
| Josephson effect for junctions which include triplet p-wave superconductors. 

One of the most significant difference of unconventional superconductors from conventional ones is the existence of 
the internal phase in the pair potential. The important role of the internal phase on the current-phase relation of the 
Josephson current has been pointed out theoretically Jl^| . In a two-dimensional junction configuration with cf-wave 
superconductors, for a fixed external phase difference between two superconductors, the component of the Josephson 
t-H ! current becomes either positive or negative depending on the injection angle of the quasiparticle fl3[| . In addition to 
^ • this, recent studies have found the formation of the zero-energy states (ZES) near the surfaces and boundaries under 
the influence of the internal phase [[l4] |l6| . The critical Josephson current Jc{T) via the ZES |J|,[D3@ is shown 
to be proportional to T — 1 . The interplay of above two effects induces a non-monotonous temperature dependence 
of the Josephson current depending on the orientation of the junctions. As for the current-phase relation, it can be 
^ decomposed into the infinite series of svn(rap) with integer n, 
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J((fi) =^2j n sm(nip), (1) 



where ip is the difference of the external phase between two superconductors 19 1. It was revealed in previous papers 
pfj|-p2| that Ji vanishes in singlet superconductor / triplet superconductor junctions independent of the orientation of 
the junctions if the spin orbit coupling is absent. This peculiar property originates from the difference in the parities 
between the two superconductors. 

Recent experiment by Jin reported a non-monotonous temperature dependence of the critical Josephson current 
£j , in Pb/Sr2Ru04/Pb junctions p3|. Since this behavior cannot be understood in terms of conventional theory of the 
Josephson effect, it is interesting to clarify the origin in terms of the influence of the p-wave pair potentials on the 
.£h ! Josephson current lying between two s-wave superconductors. In this paper, we present a theory of the Josephson 
current in a s-wave/ p-wave/s-wave (s/p/s) trilayer superconductor junction where the transition temperature of 
p-wave superconductor T p is lower than that of s-wave superconductor T s . 

In the following, we assume an s/p/s junction in the clean limit with semi-infinite s-wave superconductor region 
(x < and x > L) where the thickness of the p-wave superconductor is L. The flat interfaces are perpendicular to the 
x-axis, and are located at x = and x — L. The barrier potential at the interface is ignored for simplicity. We also 
assume that the Fermi wave number Uf and the effective mass m are equal in the all three regions. To express the 
nearly two-dimensional Fermi surface, the z component of the Fermi momentum Ufz is restricted to the region given 
by —8 < sin -1 {kp z /kp) < S. The two-component wave function ty(x) of the quasiparticles is given as the solution of 
the Bogoliubov-de Gennes (BdG) equation following the quasi-classical approximations p4j . In this approximation, 
the effective pair potentials felt by the quasiparticles depend only on the direction of the motion of the quasiparticles. 
The pair potential of the system is assumed to be the stepwise form pH] (see Fig.l) 
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In the above, 9 is the polar angle and <fi is the azimuthal angle in the x-y plane. For the p-wave pair potential, we 
assume an unitary state where only A^i(x,9, cf>) [= A^(x,9, </>)] is non-zero for simplicity [ pp[ . The effective pair 
potential felt by the quasiparticles changes its sign depending on the spin index in the s-wave region. While in the 
p-wave region, the pair potential has the same sign regardless of the spins of quasiparticles. The external phase 
difference between the p-wave superconductor and the s-wave superconductor at each interfaces is ip/2 as to satisfy 
the current conservation law pj| . 

Suppose an electron-like quasiparticle (ELQ) is injected from the left superconductor with the injection angles 9 and 
4>. In the present case, the possible reflection process is only Andreev reflection because of the absence of the barrier 
potential at the interfaces. The coefficients of the Andreev reflection is determined by solving the BdG equation under 
the following boundary conditions: 
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The obtained coefficients depend on the direction of the spin of the injected ELQ. For the calculation of Josephson 
current, we will extend the formula by Furusaki and Tsukada for s-wave superconductors. In this formula, the 
Josephson current is expressed by the generalized Andreev coefficients ai(2)(y, 9, <j)) which are obtained by the analytic 
continuation from E to iui n in a^) (<p, 9, <f>) where E is the energy of the quasiparticles measured from the Fermi energy 
Ep and w n = 27rfcsT(n + 1/2) is the Matsubara frequency. Then, the Josephson current J((p) is given as 
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B ± = e^T ns T np (l - e~ 2XW ) ± (1 + T 2 e~ 2XW ), 
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Here Rn denotes the normal resistance, with X = L/^ p and W = Sl np /{\ A p (0) | sin 9 cos cj>) . The coherence 
length in the triplet superconductor £ p is given as £ p = hkp/ \ mA p (0) |. In the present paper, we select two 
kinds of unitary p-wave pair potentials belonging to the two dimensional E M representation which are expressed as 
A U (T, <?,</>) = A p (T,e,(p) = Ap(T)sin6l(sin0 + cos</>) and A p (T)e^sm(9. Hereafter, we will call these as E„(U1) and 
E U (U2) states, respectively. The quantity 7 which is expressed by the internal phase of the triplet superconductor is 
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given by 7 = 1 for E U (U1), and 7 = e 1 ^ for E U (U2). By substituting these pair potentials in the above formulas, we 
can obtain the Josephson current J{<f) and Jc(T). 

In the following calculation, the transition temperatures T p and T s are set to be 1.0K and 8. IK, respectively, in order 
to simulate Pb/S^RuC^/Pb junctions. To express the two-dimensional features of the Fermi surface, the quantity 5 
is chosen as 0.17T. The temperature dependence of A S (T) and A P (T) are assumed to obey the BCS relation. Figures 2 
and 3 show the calculated Jc(T)'s for E„(U1) and E„(U2), respectively. With the decrease of the T, the characteristic 
temperature where Jc(T) begins to decrease is given by T p . With the further decrease of T, Jc(T) increases again. 
Apparently these features are quite distinct from those for conventional s-wave superconductors. To understand this 
behavior intuitively, it is useful to present analytical forms of K(ip, 9, <j)) in the two limiting cases. For T > T p , where 
A p (T) = 0, is satisfied, K{if, 9, 0) is reduced to be 

°> 0) = E "» 20: 1 I FLe-^e- + 1 p / ' (12) 
which reproduces that of in an ordinary SNS junction [Q. While for the case L — > 00, K((p, 9, (f>) is calculated as 

1 f -27*r Tts (i + rL)r^sin^ 

(i - r 2 nXL ) 2 + 4 7 * 2 rLr^ sin%/2) 



This result coincides with that of the s-wave / p-wave junction [ pp[ . It should be emphasized that a contribution to 
J {if) from the former (referred to as s/s-coupling) is positive and the latter is negative (referred to as s/p-coupling) for 
< if < ir. This means that the junction changes to so-called "7r-junction" when the s/p-coupling becomes dominant 
compared to the s/s-coupling. Since these two coupling terms have different T and L dependence, the competition 
of these two determines the behavior of the junction. 

Next, we will examine the thickness (L) dependence of the Josephson current. For sufficiently small L (L -C £ p ), J{f) 
remains positive independent of temperatures because of the strong s/s-coupling (see upper panel of Fig. 4). When 
L ~ £p> the junction is also governed by the s/s-coupling for T > T p . However, with the decrease of temperatures, the 
contribution from s/p-coupling term [Eq. (|13|)1 is rapidly enhanced. The sign changes in J {if) occurs when the s/p- 
coupling becomes dominant at sufficiently low temperature. At temperature just below T pi a crossover occurs between 
these two terms. This is the origin of the non-monotonous temperature dependence of Jc(T). While for L ^> £ p , 
the s/s-coupling at T > T p is weakened because it is governed by the factor exp[— LT/{%vp)\. As the temperature 
is lowered below T p , the junction easily transits to the 7r-junction Then the region, where dJc(T)/dT > is 

satisfied, becomes narrower with the increase of L. Thus the non-monotonous temperature dependence is reduced 
and becomes undetectable. 

In this paper, we present a theory of the Josephson current in a s/p/s junction and obtain the non-monotonous 
temperature dependence of Jc(T). This effect originates from the competition between the positive contribution of 
K(tp,9,<p) to J {if) due to the conventional s/s-coupling and the negative one due to the s/p-coupling. The present 
effect is quite distinct from the non-monotonous Jc{T) predicted in d-wave superconductor junction in the 

point that the present effect is insensitive to the orientation of the junction. Although we have assumed unitary state 
for Sr2Ru04, it is also possible to perform the similar calculations for several non- unitary pair potentials ]q,p|,p7t . 
However, only the monotonous increase of Jc{T) is expected because of the absence of the Josephson coupling between 
s-wave and these non-unitary pair potentials. It is important to note that the present result is consistent with the 
recent measurements of Josephson current in Pb/Sr2Ru04/Pb junction ]23]| . We believe that the consistency strongly 
supports the unitary triplet p-wave paring states in Sr2Ru04. Since the difference of the parity plays an essential 
role in this effect, the present trilayer junction configuration can be applied to a novel phase-sensitive test for the 
unconventional pairing states. For example, if E u i state |2g] is also realized in UPt3, similar effect will be detected in 
Pb/UPts/Pb junctions. Throughout this paper, we ignored the barrier potential at the interface for the simplicity. 
However, the essence of the physics will not be changed even if we will take into account of this effect. Details will be 
presented in the forthcoming paper. 
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FIG. 1. Schematic illustration of s-wave superconductor / p-wave superconductor / s-wave superconductor (s/p/s) junction. 
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FIG. 2. Normalized critical Josephson current as a function of normalized temperature for E U (U1) state. T p /T 3 — 0.13. X 
is 0.5 and 1 as indicated in the figure. 



FIG. 3. Normalized critical Josephson current as a function of normalized temperature for _E„(U2) state. T p /T s — 0.13. X 
is 0.5 and 1 as indicated in the figure. 
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FIG. 4. Normalized Josephson current as a function of the phase difference ip for E U (XJ1) state. T p /T s — 0.13. X is 0.5 
and 1 as indicated in the figure. For X = 0.5, a:T/T a = 0.15, h:T/T s = 0.13, c:T/T s = 0.11, and d:T/T a = 0.03. For X = 1, 
a:T/T s = 0.15, b:T/r s = 0.13, c:T/T s = 0.12, and d:T/T s = 0.09. 
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